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Using an exact solution with the hyperbolic-tangent profile for the electron density in the sun, 
which is developed by Notzold and later modified by Nakagawa, we have analyzed the solar neutrino 
problem. An interference term in their approach is correctly taken into account. Combining the 
hyperbolic-tangent profile with the BP2000, we obtain a phenomenological fitting in the analytic 
form. Combining recent observed results for survival probability Pfve —> Vn) by the SNO, SK, 

SAGE, Gallex, GNO and Homestake Collaborations, we obtain a large mixing angle (LMA) whose 
figure is looking like a shoulder. 

PACS numbers: 26.65.-|-t, 14.60.Pq 


I. INTRODUCTION 

One of the most interesting subjects in elementary par¬ 
ticle physics is the solar neutrino problem. As a possible 
explanation for it, the MSW mechanism seems the best 
solution iiH- Very recently the SNO Collaboration 
has reported that the survival probability of Ve from the 
sun is P{ve —> Ve) = 0.348 ± 0.03, and the large mix¬ 
ing angle (LMA), Am^ = — m\ = 0.35 x 10“^(eV^) 

and 6 = 32°, is favor for the explanation of the global 
measurements!^, |], ^ (mi, m 2 are masses of the neu¬ 
trino mass eigenstate vi, V 2 , respectively. Hereafter 
m 2 > mi is assumed. 9 is the mixing angle defined by 
Ue = cos 9 vi -\- sin 9 V 2 and = — sin 9 vi -\- cos 9 V 2 . )■ 

On the other hand in 1988 Notzold proposed an ex¬ 
act solution for the solar neutrino oscillation, using the 
hyperbolic-tangent profile for the density of electrons de¬ 
fined by Ne in the sun|^ : 


reduced about 26% in the figure of /E^ vs. sin^ 9, as 
Eq.(||) is used for the electron density in the sun[|j. His 
result with the standard solar model (SSM) in 1988 [|| 
is shown in Fig.0(b) (V = 1.718 x 10^ I = 0.155, ro = 
-0.105 for SSM in 1988). 

In this report we would like to consider the BP2000[^ 
with Eq.(|^) and the solar neutrino problem based on 
Notzold-Nakagawa’s approach including the interference 
term. Our phenomenological fitting is shown in Fig.|^ 
{V = 1.740 X 10“^, I = 0.148, ro = -0.115 for BP2000). 

In 2nd section, their theoretical formulas a la Nak¬ 
agawa are explained. The explicit interference term is 
shown. In 3rd section, we estimate the allowed regions 
(Am^ and mixing angle) for survival probabilities ob¬ 
served by the SNO, SK, SAGE, GNO, Gallex and Home- 
stake Gollaborations[Q ||, ||, |^, |^, In 4th 

section, concluding remarks are presented. 


Ne = Ao[I - tanhy]. 


II. THEORETICAL FORMULAS WITH THE 
HYPERBOLIC-TANGENT FOR Ne 


where Nq = V / {RsV^Gf), Rs and Gf denote the ra¬ 
dius of sun Rs = 6.96 x 10®km and the Fermi coupling 
constant Gf = 8.917 x 10“®(GeV • fm^). Moreover V, r 
and I denote the magnitude of the density at the center of 
the sun, r = R/Rs and I = l/Rs, respectively. Notzold 
has assumed r —>■ r -i-rg = r — 00 , a priori (See Fig.|^(a)). 

Later, Nakagawa has proposed a modified expression 
for Nf. as 


Ne = Ao[l - tanh — 


( 2 ) 


To explain theoretical formulas first derived by Notzold 
and later modified by Nakagawa, we briefly describe their 
framework. The coupled Schrodinger equations for the 
solar neutrino oscillation are given as 



/ m? cos^ 0 -I- mn sin^ 9 , \ 

( ^- +V2GFNe]iye 


Am? 
' 2E 


sin6*cos6*r'^, (3) 


where ro is an adjusted parameter which reproduces the 
Ne in better way than Eq. (0)- Nakagawa has stressed 
that the result by Notzold, the magnitude of Am?/E^, is 
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. dv,. 


Am? m? sin^ 9 -\- m? cos^ 9 

+- - - 


( 4 ) 


Introducing the following notations, ct = Rsr, Mg = 
Rs{tTi\ m\) /{4,E), RsAm?/{2E) = x/B = xq with 
B = i?(MeV)/Am^(eV^) and x = 1.764 x 10®, we use 
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FIG. 2: The electron number density from BP2000, which is 
fitted by Eq.(^). 


we have the following differential equation as 


d'^v 


ilvUiry + {^VU{r)-Gr + H^ 


v = 0, (8) 


where U{r) = 1 —tanh[(r + ro)/?]. Changing the variable 
r to y and assuming the factorized form for v(r), 


y = 


1 + exp 



(9) 


[/(r) = 2(1 - y), (10) 


FIG. 1: (a) The electron number density from SSM1988 in 
Ref.|^, which is htted by Em®. See Fig.2 of Ref.0. (b) Ns 

is htted by Eq.(||). See Ref.||]. y{r) = (1 - yYy''f{y), (11) 


the following abbreviations 

xo cos 29 


G = 


H = 


xo sin 29 


Moreover, exchanging the variable from i/g to v(r), 

Ve{r) = v{r)expiy-i J [Mq + ^VU{r)]dr^ , 


(5) 


( 6 ) 


(7) 


we obtain the following hyper-geometric equation 

+ (.0. + ^^ + ^)y}^ - = o, (12) 

ay ay 


where 


Q — jj, ly 


b — fi u \ — A, 


(13) 


(14) 
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c = 2iy + 1, 


fi = i^-VG^+H\ 


(15) 


(16) 


(17) 


of the oscillation. We have obtained following formula, 
after integration, 


-Ks 


10.16 X lO^ALAp 

. A.08 X 10^ ^ ^ , 

X sin(- ALAp) cos 

Rs 


VPc(1 - Pc) cos 29m sin 29 
5.08 X 10^ 


Rs 


-{L 2 + Li)Ap 


(27) 


where Li and L 2 are the minimum distance from the sun 
to earth and the maximum one, respectively, and 

AL = L2-Li. (28) 


A = *^F. (18) 

The general solution of Eq. dil) is expressed by two in¬ 
dependent hyper-geometrical functions. F(a, 6;c;y)_ is 
expressed by replacing v —v in o, 6, c of F(a^ b; c; j/)+. 

fiy) = CiF{a, b\ c; y)+ + C 2 y~‘^''F{a, b\ c; y)-, (19) 


III. ANALYSIS OF SOLAR NEUTRINO 
PROBLEM WITH EMPIRICAL VALUES BY 
MEANS OF EQ.(20) 


The data on the solar neutrino, the survival probability 
P{ve —> Ve)^ reported by the SNO, Gallex, GNO, SAGE, 
Homestake Collaborations are shown in Table |. 


where Ci and C 2 are the integral constants. Using an or¬ 
dinary procedure with the initial condition, Ve{y = yo) = 
1 and v^{y = ?/o) = 0, we obtain the following formula 
for the survival probability 

= Pi cos^ 0 -I- (1 - Pi) sin^ 9 
— \/^c(l — Pc)cos29m sin20cos(2.54— —L + S). (20) 


TABLE I: Data on the solar neutrino|^, ^ |^, |^. {E^} is 
the average neutrino energy. 


Exp 

(E? 

P{Ue Ve) 

SNO 

~ 8 MeV 

0.348 ± 0.029 

SAGE 

~ 0.8 MeV 

0.54 ±0.06 

Gallex -t GNO 

~ 0.8 MeV 

0.56 ±0.07 

Homestake 

~ 8 MeV 

0.34 ±0.03 


where 


Pi = Pc sin^ 9m + ii- Pc) cos^ 


( 21 ) 


Pr = 


cosh(7r^U) — cosh[7r/(Ap — Aq)] 
cosh[7rl(Ap -I- Aq)] — cosh[7rZ(Ap — Aq )]' 


( 22 ) 


cos 20m = 


G-U(l-yo) 


^[G-V{l-yo)]^ + H^’ 


(23) 


Ap = ^Am\ 


(24) 


Using empirical values in Table || and Eq.(^) with 
(Py) = 8MeV and (E,^) = O.SMeV, we obtain two 
allowed regions (AR) which are shown in Fig.||(a). 
The small mixing angle (SMA) region is located in 
(tan^ 9, Aw?') = (0.8 x 10“^, 10“^) and the large mixing 
angle (LMA) region is located in (0.5, 10“"^ ~ 4 x 10“^). 
Combining the global measurement (including the day- 
night effect by the SNO Collaboration and by the SK 
Collaboration) and our result, we obtain Fig.||(b). The 
shoulder-like region nearly coincides with the large mix¬ 
ing angle (LMA) with C.L.99.73%. (See Ref.Q) 

As is seen in Fig.| our results based on Eq.(||) and 
Eq.(p0[) is partially consistent with the LMA solutions. 


IV. CONCLUDING REMARKS 


Aq = a/U 2 _ 2VApcos29 + Ap^, (25) 


yo = - — ■ (26) 

1 -h exp(-2y) 

In Refs.|0 dl the third term of the right hand side, the in¬ 
terference term, in Eq.(pO|) is assumed to be zero because 


1. As is seen in Fig.| Eq.(||) is useful to reproduce 
the electron density based on BP2000. 

2. Using Eq.(d^, we have examined the AR expressed 
by (Am^, tan^0). In the present calculation we 
have added the third term Eq.(p^) to Notzold and 
Nakagawa’s formula. As compared with result of 
Refs.Qdj, we have additional AR: The contribu¬ 
tion is seen in the most right curve and the bot¬ 
tom region with oscillations near (Am^, tan^ 9) = 
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(0.2 ~ 5, 1 X 10“^'^ ~ 1 X 10“^^).The main reason is 
attributed to Eq.(p^. Moreover, we can compare 
our results with the exponential profile, for exam¬ 
ple, see Fig.2 of Ref. M ■ A similar allowed region 
(0.1 < tan^ 9 < 1.0) is observed. 


3. From the empirical results in Table g and Fq.(20), 
we have obtained two AR’s, the LMA and SMA. 
Combining the AR’s reported by the SNO, the SK 
Collaborations and our Fig.^(a), we have obtained 
Fig.^(b). This is fairly well consistent with the 
LMA reported by the SNO Collaboration. 


4. Moreover, we have to add the following fact; For 


the figure of (sin^ 0, Am^/Ai,) with the BP2000, 
the result of calculation a la Ref.0, the magni¬ 
tude of is about 20% larger than that 

of present calculation using values reported by the 
SNO Collaboration^. 
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FIG. 3: (a) Am^ vs. tan^ 6. Long dashed and solid 

lines are obtained from empirical values by SNO Collab¬ 
oration {P{ve —> I'e) = 0.348 db 0.03), and those by Ga 
experiment(P(i/e ^ r'e) = 0.55 ± 0.06) and assuming {Ev) = 
8MeV and {Ev) = 0.8MeV, respectively, (b) Enlarged figure 
of (a). The dashed circle with C.L. = 90% and dotted circle 
with C.L. = 99.73% are given in Ref.^. The SMA is excluded 
by observations by the SNO and SK Collaboration. AR is an 
abbreviation for allowed region. 






















